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Abstract 

The dynamics of a physical gel, namely the Low Molecular Mass Organic Gelator Methyl- 
4,6-O-benzylidene-a -D-mannopyranoside (a-manno) in water and toluene are probed by neu¬ 
tron scattering. Using high gelator concentrations, we were able to determine, on a timescale 
from a few ps to 1 ns, the number of solvent molecules that are immobilised by the rigid net¬ 
work formed by the gelators. We found that only few toluene molecules per gelator participate 
to the network which is formed by hydrogen bonding between the gelators’ sugar moieties. 

In water, however, the interactions leading to the gel formations are weaker, involving dipolar, 
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hydrophobic or n — n interactions and hydrogen bonds are formed between the gelators and the 
surrounding water. Therefore, around 10 to 14 water molecules per gelator are immobilised 
by the presence of the network. This study shows that neutron scattering can give valuable 
information about the behaviour of solvent confined in a molecular gel. 


Introduction 

Physicals gels formed by low molecular mass organic gelators (LMMOG) represent very interest¬ 
ing materials, both for fundamental physics as well as applied science. They are composed of a 
rigid network formed by the gelators, in which is trapped a large quantity of solvent. Gathering 
properties of both liquids and solids, they find applications in diverse domains, such as temperature 
sensors or drug delivery- -. They present a high potential for functionalised nanomaterials, thanks 
to the versatility of functionalized molecules that can be used as gelators- - -, and the supramolec- 
ular organisation that can be tuned by changing the solvent. 

The same gelator can indeed present various structural organisations, like ribbons, vesicles or 
even nanotubes 9 ’ 10 ', that may arrange in periodic uni- or bidimensional structure. The particular 
arrangement adopted is defined by a subtle interplay between different forces like hydrogen bonds, 
hydrophobic forces, n — n interactions or dipolar forces, so it will depends on the solvent and will 
be often restricted to a particular temperature range. The gelation of new molecules, even after in- 
tensives research based on steric effects and crystalline structure remains poorly predictable 81112 . 
Moreover, applications suffer several limitations, due to the metastability of the gels that may col¬ 
lapse after some delay, spanning from the hour to the month or longer. Although having certain 
organisation, it is known that these systems are rather disordered and can show several relaxation 
modes related both to the gelators and the solvent dynamics. 

If structural information on LMMOG obtained by Transmission Electron Microscopy (TEM) 
or small angle scattering is often available, this is not the case concerning the dynamics. Their 
stability can be evaluated from the gelation temperature and melting enthalpy— - - 5 and compared 
to solvent parameters characterizing its polarisability or hydrogen bond ability 14 - 1617 , leading to 
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a phenomenological description of the gelator-solvent interactions. Infra-Red spectroscopy can 
be used in order to investigate the strength of hydrogen bonding 1 -, related to the stability of the 
network. The dynamics of both gelators and solvent molecules are important contributions to the 
gel stability, as a source of entropy. The solvent dynamics should be influenced by the presence 
of the gel, specially if some solvent molecules participate into the rigid network. In that case, 
several solvent populations may be distinguished. The gelators mobility is also greatly affected 
by the gel formation. Mainly NMR-based techniques have been used to quantify the mobility of 
the gelators and the solvent in both the sol and gel phases^ 1 , on timescales slower than tens 
of nanoseconds. The solubility of the gelators in the solvent, which is small but finite, leads to 
a dynamical exchange of gelators between the rigid network and the liquid part. The network 
is therefore a dynamical assembly—. The life time of a gelator in the rigid network is expected 
anyway to be longer than the nanosecond—’ 22 , so gelators appear immobile on a shorter timescale. 

Given the typical gelators concentration and sample slight turbidity, the characteristic size of 
the pores is estimated to be larger than a few hundred nanometres. A large part of the solvent is 
therefore expected to present a bulk-like behaviour. However, the interaction between the gelators 
and the solvent molecules leads to an organisation of the solvent close to the rigid parts, as well as 
a slowing down of its dynamics 23-25 . Such behaviour is commonly observed for liquids confined 
in matrices. Water confined in porous silicas is possibly the best characterized example of this 
kind of system and it has been shown that even in pores of ~ 2 nm of diameters, only 2 water 
layers are influenced by the matrix, and the third layer presents already a bulk-like behaviour 26 . 
A similar behaviour is found for confined solvents in hard or soft confinement matrices—. The 
dynamics of solvent is expected to exhibit a continuous transition 27 between the most hindered 
(interacting) solvent molecules and the bulk-like ones. In practice, this behaviour is experimentally 
characterized by the distinction of several populations, varying between 2 and 4 depending on the 
experimental techniques 2 -. In molecular gels, the mobility of the solvent has been investigated by 
NMR 18-22 on timescales slower than 20 ns. Quasi-elastic Neutron Scattering operate on a faster 
timescale and is therefore particularly well suited for characterising the faster solvent dynamics, 
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but it has barely been used to investigate the dynamics in gels, as the solvent contribution usually 
overwhelms the whole signal 29 . 

In this context, we undertook the study of the microscopic dynamics in a model physical molec¬ 
ular gel, with the goal of characterising the molecular interactions and understanding the formation 
and stability of the gel phase. As model system, we chose to investigate a gelator that belongs 
to a sugar-based gelator family first reported by the group of Shinkai and collaborators, methyl- 
4, 6-O-benzylidene-a-D-monopyranoside 30 . It is one of the simplest organic molecular gelators 
family described in the literature, made of a saccharide carrying an aromatic benzylidene group as 
schematized in the Schemed] The methyl-4,6-O-benzylidene-a-D-glucopyranoside (a-gluco) and 
methyl-4,6-O-benzylidene-a-D-mannopyranoside (a-manno) derivatives present the best gelation 
properties. They only differ by the orientation of one of the OH groups (see Schemed]), but the 
a-manno gelifies protic solvent (water) as well as apolar organic solvents (toluene), while the a- 
gluco does not gelify water, a-manno is therefore a small molecule carrying not more than the 
necessary functions for an amphiphile gelator. Its compact shape enables to use simple models for 
its dynamical relaxations (rotational and translational diffusions) and possible interactions with the 
solvent. 

Like most gelators, a-manno adopts different supramolecular conformations depending on 
the solvent, forming more or less complex assemblies. For example, a remarkable structure is 
the one formed in p-xylene, in which fibres of a-manno are organized in an hexagonal structure 
with characteristic distances of 60 A, involving solvent molecules in the rigid network- 1 . As a 
general trend, the gelators assemble in water through the stacking of the benzylidene parts and 
form hydrogen bonds with the solvent while in apolar or weakly polar solvents, the gel formation 
is driven by intermolecular hydrogen bonds between the gelators and the aromatic groups are 
solvated by the liquid part 10 - 32 . 

With the goal of extending the investigated timescale of the dynamics in molecular gels, we 
performed quasi elastic neutron scattering on the model system a-manno in water and toluene. 
We aim at characterizing how the solvent dynamics are affected by the gelification and how it 
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compares in different kind of solvents. Combining Time Of Flight and Backscattering data, we 
are able to span the dynamics over four orders of magnitude, from tenth of picoseconds up to 
the nanosecond in order to characterize the various solvent populations. In all our measurements, 
given the characteristic timescales previously discussed, the gelators are therefore considered to 
be immobile for translation. After the overview of the Material and Methods, we will present the 
experimental results and analysis, and discuss them in regard of results from the literature obtained 
by NMR techniques. 

1 2 

Scheme 1: methyl-4,6-O-benzylidene-a-D-glucopyranoside (1) and methyl-4,6-O-benzylidene- 
a-D-mannopyranoside (2) 




Materials and Methods 

Sample preparation, a-manno was purchased from Orgentis Chemicals and dried 

under vacuum. Organic solvents were bough from Sigma-Aldrich, distilled at 80°C and dried 
over molecular sieve, a-manno powder was mixed with water or solvent and heated to 80°C, 
and the gels were formed during cooling at room temperature. The Table Q] describes the different 
samples used on both instruments. In order to check the concentration dependence of the net¬ 
work constituents, we performed the measurements on IN5 (see next paragraph) on three different 
contributions of a-manno toluene gels. In all samples, the residual solubility of the gelators was 
always lower than ~ 10 1 g / L, and therefore neglected as a contribution in the liquid part. We 
checked before and after all measurements that the gels did not collapse. 

Neutron Scattering. Quasi elastic neutron scattering was performed at the Institut Laue- 
Langevin on the spectrometers IN5 33 and IN16— on fully protonated samples. On IN5, the quasi 
elastic signals were measured using a wavelength of A =5 A, providing a resolution of 100 /leV 
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Table 1: a-man no gels used in the QENS experiments and contribution of each component to the 
total neutron scattered intensity. 


instr. 

solvent 

[c] g/L 

numb, solvent/gelator 

I gelators 

I solvent 

IN5 

toluene - T1 

51.2 

52.1 

0.04 

0.96 


toluene - T2 

103.6 

25.7 

0.08 

0.92 


toluene - T3 

145.0 

18.4 

0.11 

0.89 


H 2 0- 1 

88.2 

177.5 

0.05 

0.95 

IN16 

toluene - T4 

139.3 

19.1 

0.10 

0.90 


H 2 0 - 2 

208.4 

75.2 

0.11 

0.89 


and the Q range [0.2-2.2] A -1 as well as a wavelength of A=10 A, providing a resolution of 10 
jUeV and the Q range [0.1-1.1] A -1 . On IN16, the instrumental resolution was 1/ieV and the Q 
range [0.2-1.9] A -1 . The quasi-elastic spectra were measured between -10 and 10 /reV, while the 
so-called fixed window scans enabled the monitoring of the (elastic) intensity integrated over the 
instrumental resolution in a temperature range of [50-350 K]. In these last experiments, the samples 
were first heated up to above the sol-gel transition, and then cooled down to low temperature, to 
avoid an eventual damage of the gel at low temperature. The sol-gel transition is indeed reversible, 
but the effect of the solvent crystallisation in the network may cause damages in the gel structure. 

The Table Q] also reports the fraction of neutron scattered intensity of both constituent of each 
sample. The calculations are based on the incoherent cross section, widely dominating the whole 
signal: G a - marmo =l43SA, <j to i uene =639.3 and o water =159.& bams, respectively for the a-manno, 
the toluene and the water molecule. The data reduction and basic corrections were performed 
using the standard procedures in the LAMP program^. The IN5 data were rebinned at constant Q 
values, but for the analysis of the elastic intensities the constant 20 binning was used, enabling the 
use of a larger dynamical range. The sample thicknesses were of 0.2 mm, ensuring a transmission 
higher than 0.87 and enabling to neglect multiple scattering. 


Experimental results 


Fast solvent dynamics. We first investigated the fast dynamics using Time of Flight measurements 
(IN5, ILL) with two different resolutions, giving access to the timescale 0.1-100 picoseconds. The 
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bulk diffusion coefficients of water and toluene at 293 K are respectively around 2.0-10“ 5 cm 2 /s^ 
and 2.7-10 5 cm 2 /s 37 . The investigated dynamical window is therefore suitable to observe such 
translational dynamics of bulk solvents, as well as dynamics 10 to 15 times slower. 



-3 - 2-10 1 -3 - 2-10 1 


energy transfer (meV) 

Figure 1: Picosecond dynamics measured by QENS on the Time of Flight IN5 at 290 K, Q= 
1.3 A -1 in the gel of a-manno in toluene (T3) and water (HiO-l), and pure toluene and water for 
comparison (bottom). The lines show the fitting to the experimental points using the equation (??) 
(red line) and the different elastic and quasielastic contributions (green lines)(see text). 


Both gels and pure solvents were measured in the same conditions for comparison. Charac¬ 
teristics quasi-elastic neutron scattering spectra of the gels and corresponding pure solvents are 
shown in Figure |T] The main difference observed between the gels and pure solvents is an elastic 
contribution that emerges in the gel as a narrow peak on top of the broad quasi elastic contributions 
arising from the liquid fraction, evidencing the presence of the rigid network. 

Both set of data, in water and toluene, were fitted with a phenomenological model described 
by the equation ?? and made of two lorentzian contributions and an elastic peak. This model is 
not the most sophisticated one that can be used to describe the dynamics of water and get a proper 
separation of the rotational and translational dynamics. However, it better suits our purpose, which 
is the differentiation of water dynamics in bulk and in gel and the comparison of the translational 
diffusion coefficient, enabling us to use the same model for all the samples investigated here. This 
theoretical form was further convoluted with the instrumental function, and no constraints were 
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applied to the parameters. 


5(2, co) = (A 0 (g) • 5(©) + Ai (G)Li(2,©) +A 2 (2)^2(2, ©)) ® *(2, ©) (1) 

All the spectra exhibit a broad lorentzian having a width of several meV with a mild 2-dependence 
that is associated to rotational dynamics, and a narrower lorentzian (0.1-1 meV FWHM) showing 
the characteristic 2-dependence of a translational diffusion process. The agreement with experi¬ 
mental data was satisfactory, and no additional contribution of intermediate relaxation time could 
improve the fit. 

Translational component. The width of the translational contribution over the whole Q range 
is represented in Figure |2]for the different samples and the two experimental configurations. In the 
three toluene samples, no difference can be observed between the gel and pure solvent, meaning 
that no slowing down of the solvent dynamics could be detected, whatever the gelator concentration 
is. More sophisticated fits using rotation-translation models, or taking extensively account of the 
solvent in the gel samples lead to the same result: no intermediate dynamics is observed between 
the bulk-like behaviour of the sample and the molecules having a relaxation time slower than 100 
ps, and that are therefore detected as immobile, contributing to the elastic peak. 

The toluene dynamics follows a Fickian behaviour, with a width of the quasi-elastic compo¬ 
nent increasing linearly with Q 2 . We can however observe a deviation of this law at Q ~ 1.2A -1 , 
corresponding to a distance of 5.3 Awhich is close to the distance between the center of mass of 
the molecules. This plateau is typical of the De Gennes narrowing, that consists in a slowing down 
of the translational dynamics when the translational distance approach the equilibrium distance 
between the molecules. This effect is observed only in the coherent neutron scattering signal (8% 
of the total scattering in the case of toluene), that increases at this momentum transfer according 
to the structure factor. The toluene dynamics can however be well approximated by a simple Fick¬ 
ian diffusion, with a coefficient of 3.2-10 _5 cm 2 .,s _1 . The fit of this coefficient is represented in 
Figure |2] (bottom). This slight overestimation with respect to values reported in the literature is as- 
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signed to the use of separated components between rotation and translation, instead of performing 
a convolution between the two components as would be required by a more accurate treatment. 

In water, the translational dynamics is modelled according to the jump diffusion behaviour—. 
The respective coefficients for pure water and hydrogels are D=2. 8±0.T10 5 cm 2 .s 1 and D=2.5 
± 0.1T0 _5 cm 2 .5 _1 , and the residence times T = 0.75 ± 0.03 ps, and z = 0.77 ± 0.04 ps. The 
translational dynamics of water is therefore sightly slower in the gel than in pure water. We note 
that the single-lorentzian character of the signal proves that the whole water is slowed down and 
not only a fraction of it that would be averaged with the other bulk-like fraction. The effect, as will 
be discussed in more details in the Discussion part, is however very small compared to what ca be 


observed on other systems. 



0,0 0,5 1,0 1,5 2,0 

Q (A 1 ) 


Figure 2: FWHM of the translational component of the solvent signal in the pure solvents and 
associated gel measured on IN5 with a wavelength of 5 (squares) and 10 (circles). Top : a manno- 
HiO (black symbols) and pure HiO (red symbols), bottom : a manno-Toluene (black symbols) 
and pure Toluene (red symbols). The fit of the jump diffusion in water samples and Fick diffusion 
for Toluene samples are represented as solid lines for comparison. The diffusion coefficients are 
discussed in the text. 


Rotational component. The large component is typical of rotation or localized motions, ex¬ 
hibiting a rather flat width and increasing amplitude when Q increases. In pure toluene, the com- 
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ponent has a full width of 2.7 ± 0.3 meV, while the corresponding component in the gel, having 
the same characteristics at small Q, tends to increase for Q larger than 1.4 A -1 up to ~ 6 meV. The 
case of water gels is very similar: while the pure water exhibits a rotational component of constant 
FWHM =1.5 ±0.2meV, the quasi elastic component width in the gel increases up to ~ 6 meV. 
The increase of the component at high Q is interpreted as a quasielastic component arising from a 
localized motion of the gelators within the network, like flip-flop motions of the sugar moieties or 
rotation of the aromatic groups. 

Elastic component. Although no influence of the network can be detected on the solvent 
mobility, quantitative information about the network composition can be extracted from the elastic 
peak, which arises from all the atoms in that system that do not diffuse or appear as immobile 
on the timescale given by the instrument resolution, i.e. t <13 ps. In order to evaluate the elastic 
contribution, the fitting of the data measured in the gels was performed using the pure solvent 
as a (scaled) contribution. This way, only an additional elastic contribution (delta function) was 
necessary to fit the data. We should note here that, as stated in the discussion of the rotational 
component, a quasi-elastic signal arising from the network should be detected at high Q. However, 
it is not possible to reliably distinguish a new rotational component in the data. The small quasi¬ 
elastic signal of the gelators is therefore drown in the solvent rotational contribution and the error 
bars. 

The total scattering function may then be written as a sum of liquid and solid contribution. As 
previously shown by the analysis of the translational diffusion coefficients, the liquid parts exhibit 
diffusion coefficients which are highly similar to those of the bulk solvent, enabling us to identify 
the liquid contribution with the pure solvent one: S(Q, (0 ) = CcS so nd(Q, (o) + (1 — <x)S so i vent (Q, co). 
The solvent part is therefore purely quasi-elastic while the solid one is the sum of elastic and 
quasielastic terms due to localised motions of gelators or solvent molecules trapped in the rigid 
network. Moreover, the total scattered intensity is conserved such that I e i + I qens = 1. The quan¬ 
tity R(Q) = ( - el,solid -| rcca ]] s the Elastic Incoherent Structure Factor (EISF) usu- 

ally defined in solids 39 . In the limit of Q —> 0, all the atoms belonging to the solid part, and 
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only them, contribute to the elastic component, as their quasi-elastic contribution drops to 0. We 
therefore get an evaluation of the solid part from the ratio R(0) — 

Gamanno “b Msolvent,immo ' &solvent \ 


|, where o sohem and <7„ 


^el,solid “b 1qens,solvent J q —q 

are the scattering cross sections and 


Gamanno ~b nsolvent,total ' ®solvent J 

nsolvent,immo and n so i vent are the number of solvent molecules per gelator, respectively immobile and 
total, assuming all the gelators contributing to the network. 

We eventually get the number of immobile molecules from the relation: 


Gsolvent ‘ ^-solvent ' 


^solvent .immo 


lei,sol id 
Iqens,solvent J Q=() 


-o, 


amanno 


Gsolvent ' ( j 
\ 1 ! 


I el,sol id 


'qens,solvent 


+ 1 


Q =o 


( 2 ) 


These quantities are shown in Figure [3] for the four investigated samples. As can be observed 
for the toluene samples, a Q dependence is observed, showing that a relaxation is present even in 
the rigid network. Analyzing this Q dependence gives access to characteristic geometry and length 
scales of the undergoing dynamics in the solid. The variation of R(0) as a function of Q can be 
fitted with the structure factor given by an isotropic rotation: 


R(Q)=A 


( sin ^(3 

V r s Q 


2 

+ B 


(3) 


where A is the amplitude, B is the intensity corresponding to the atoms that do not participate to 
the relaxation, and r s is the radius of rotation. The equation ?? gives a reasonable description of 
the data and enables, beside extracting a radius of rotation, to extrapolate the data towards Q — 0 
and get the correct value for R( 0). 

The number of solvent molecules n so i ventjimmo needed to reproduce R(Q) can then be extracted 
from the data. This number is the number of solvent molecules that are included in the rigid 
assembly. All the data are collected in Table El showing that an average of 6.2 ± 2.5 toluene 
molecules and 10 ±3 water molecules are immobilized on a timescale of 13 ps. Concerning the 
different toluene gels, the number of trapped solvent molecules does not, within the error bars, 
vary significantly with the gelator concentration. We will therefore assume the network constitu- 
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tion as concentration independent, and consider only the average number of immobilized toluene 
molecules. The simple isotropic rotation model, although providing a poor fitting of the oscilla¬ 
tions at high Q due to the superposition of other relaxation processes, enables to extract a mean 
rotation radius. A value of r s = 3.8 ±0.2 A was extracted for the three samples, in good agree¬ 
ment with the rotation radius that could be expected from the size of the benzylidene moieties of 
the gelators or of the toluene molecules trapped in the network. The water sample, instead, do not 
exhibit any Q dependence of the elastic intensity, indicating that no dynamics of the solid part of 
the sample are detected on this timescale. 

Table 2: a-manno gels used in the QENS experiments, ratio between elastic and total intensity 
R( 0) measured on IN5, and number of immobilised solvent molecule per gelator on 13 ps. 


sample 

solvent/gelator 
num. mol. 

m 

norm. 

solvent immobilised 
num. mol. 

T1 

52.1 

0.20 ± 0.02 

8.8 ±3.8 

T2 

25.8 

0.24 ± 0.02 

4.5± 0.8 

T3 

18.5 

0.35 ± 0.02 

5.1± 2.6 

H 2 0- 1 

180.2 

0.1 ±0.02 

10 ±3.0 



Figure 3: Ratio between the elastic contribution and the total intensity measured (integrated over 
an energy range of [-5,1.8] meV) at A = 5 A on IN5. The fits to the toluene data represent the 
rotational EISF (see text for details). The inset shows the same data for the water sample where no 
significant Q dependence is observed. 


Sample fraction immobile up to 1.3 ns. Quasi elastic measurements on the backscattering 
spectrometer IN 16 (ILL) showed that no relaxation with characteristic time between 100 ps and 
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1.2 ns is present in the gel phase at 300 K. But the elastic intensity, recorded as a function of 
temperature in a so-called fixed window scan gives quantitative information on the solvent fraction 
that is effectively embedded in the rigid network--^, in a similar way the elastic contribution did 
for the IN5 data. In such a scan, a calibration of the intensity can be performed by measuring at low 
enough temperature, typically 50 -250 K, where the whole sample is frozen. In this temperature 
range, the whole sample contributes to the elastic intensity, which decreases linearly with the 
temperature following a slope given by a pseudo-Debye-Waller factor, i.e. decrease of the intensity 
due to thermal agitation. The intensity is therefore proportional to the total number of molecules 
present in the sample, pondered by the suitable cross sections. The extrapolation of this slope at 
300 K equals the intensity that would be measured if the whole sample was still solid, i.e. did not 
transit from the solid to the gel state (melted solvent). At 300 K, the bulk-like liquid fraction of 
the sample moves faster than the resolution and therefore do not contribute to the elastic peak any 
more. Like in the previous analysis, we assume that all the gelator molecules participate to the 
rigid network, for which we can calculate the elastic intensity. The difference of intensity between 
the experimental one and the one that would be given by the gelators can therefore be attributed to 
the solvent molecules that are immobilized together with the gelators. 

, Idw G (Xmanno T ^ solvent ' G.solvent \ , , . , ...... 

We have-= --r, where I exp is the measured intensity at 300 

lexp [Gamanno ± Insolvent,immo ' Gsolvent\ 

K and low is the intensity extrapolated at 300 K following the Debye-Waller slope. The data are 
represented in Figure |4] and the results are given in Table [3} 

Table 3: a-manno gels used in the IN16 experiments, IoW, I e xp and number of immobili z ed 
solvent molecule per gelator. 


sample 

num. solvent/gelator 

Iexp(300K) 

I D.W. 

num. solvent immobilised 

T4 

22.0 

1.25±0.25 

8.6±0.25 

1.23 ±0.5 

H 2 0 - 2 

75.2 

2.60±0.25 

9.60±0.25 

14 ±2 


At 300K, this procedure gives n so [ verit j mmo = 1.23±0.5 immobile toluene molecules in the a- 
manno/toluene sample, and n so i vent j mmo =14 ±2 in the a-manno/water gel. The result for water 
is slightly higher than the value obtained from the IN5 data corresponding to a shorter observa¬ 
tion time ( 13 ps for IN5 and 1.3 ns for IN16), which obviously is unphysical. However taking 
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into account the large error bars associated with the procedure followed to estimate the number of 
immobilized molecules, both results agree reasonable well and allow us to conclude that approxi¬ 
mately 10-14 water molecules associate strongly to the gel network and exhibit therefore solid-like 
dynamics instead of the bulk behavior found for the remaining water. In the case of toluene, the 
error bars are also large, but it seems clear that the number of immobilized molecules on a time 
scale of 13 ps decreases when the observation period is increased by two orders of magnitude. 



Temperature (K) 

Figure 4: Elastic scan as a function of temperature measured on IN 16 (instrumental resolution 
1 /leV) and integrated over the full Q range. The line indicates the extrapolation of the pseudo- 
Debye-Waller intensity toward high temperature, if the sample did not melt. 


Discussion 

The solvent dynamics measured by Quasi Elastic Neutron Scattering is relevant of the different 
nature of both gels. We indeed showed that only a few toluene molecules are trapped in the 
network formed by the a-manno even on the picosecond timescale, while 10 to 14 water molecules 
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are immobile up to a few nanoseconds in the hydrogel. The different behaviour between water and 
toluene can be rationalized in the different nature of the interaction of both kinds of solvent with the 
gelator, already responsible for the different structural organization of the gels. In the case of water, 
H-bonds are formed between H 2 O and the polar groups of a-manno. They are relatively strong and 
therefore fix water in a stable manner, blocking their dynamics so no quasielastic signal from those 
water molecules is observed and they remain attached to the gel for times of the order of several ns 
or larger. For toluene, the interaction between the solvent and the gelator is likely through n — k 
interactions between the aromatic rings of a-manno and toluene. Due to steric effects this makes 
that the number of fixed solvent molecules is smaller and that the dynamics of those molecules are 
not completely frozen, so that they retain some rotational freedom that produces the QENS signal 
assigned to the solid and characterized by the EISF shown in Figure [3} In this case, the weaker 
character of the interaction between solvent and gel makes that toluene can escape from the gel 
in a time scale of several hundreds of ps, so that the number of immobili z ed molecules is 5 on a 
time scale of about 13 ps and only 1 molecule when the observation time is 1.3 ns. In principle, 
we could expect that such escaping motion should appear as an additional relaxation slower than 
the bulk translational diffusion. However the relatively small number of molecules participating 
in such motion (3-4 per gelator vs a total number of toluene molecules between 18 and 52 per 
gelator, depending on the concentration) together with the instrumental limitations (both in terms 
of resolution and energy range, as well as in terms of statistical accuracy), make this difficult to 
observe. 

It is interesting to compare the present hydrogel with sugar solutions, in which recent studies 
of hydration water dynamics show a slowing down of the dynamics with respect to pure water. 
The retardation factor varies depending on the study and the sugar concentration, and is estimated 
to about 3-6 in mono- and disaccharides 42 ’ 43 , and does not exceed ~ 15 for the first hydration 
layer—. In our case, indeed, a slight slowing down of the water dynamics is observed by Time 
of Flight measurements, that could arise from the average of the bulk-like water and slowed down 
hydration water. However, the measured number of ~ 10 to 15 molecules with dynamics slower 
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than 1.2 ns indicates that the water molecules are more strongly involved in the rigid network than 
usual hydration water. 

The group of Tritt-Goc measured the diffusional dynamics of solvent in very similar gels 
made of methyl-4,6-0-(p-nitro-benzylidene)-a-D-glucopyranoside as gelator, either in water or 
chlorobenzene 20 - 45 . This gelator is also an amphiphilic one, but dipole-dipole interactions are 
more important that in the a-man no thanks to the nitro-substitute group. These forces act in favour 
of the gelation process, enabling gelation at lower concentration, and more stable gel phases than 
with the a-man no. In water, the translational diffusion coefficients of gels with concentration 
ranging between 1-5 % wt. is similar to the one of bulk water, over distances of a few microns. 
They could not distinguish any component characteristic of water slower than the bulk 45 . Similar 
studies, performed on a gel composed by the same gelator but in chlorobenzene, a slightly polar 
and aprotic solvent 20 , in the same concentration range, proved a significant slowing down of the 
solvent dynamics close to the gel surfaces. They were able to estimate to four the number of sol¬ 
vent molecules per gelator immobili z ed on the probed timescale, i.e. between 25 ns and 100 /is 
(10kHz and 40 MHz). Assuming both studies may be related, in sight of the slight diversity of 
gelators and solvent, we can conclude that exists an intermediate dynamical range on which the 
water molecules eventually start moving. 

In conclusion, this study demonstrates that Quasi-Elastic Neutron Scattering can be used to 
characterize the solvent dynamics in a gel, on timescales complementary to NMR. The combina¬ 
tion of both techniques indicate that several populations of solvent can be distinguished on differ¬ 
ent timescales: a few highly hindered solvent molecules with relaxation times in the microsecond 
range, and another population moving on the sub-nanosecond range. 
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